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Abstract 

The 0(n) loop model on the honeycomb lattice with mixed ordinary and 
special boundary conditions is solved exactly by means of the Bethe ansatz. 
The calculation of the dominant finite-size corrections to the eigenspectrum 
yields the mixed boundary scaling index and the geometric scaling dimensions 
describing the universal surface critical behaviour. Exact results follow in the 
limit n = for the polymer adsorption transition with a mixed adsorbing and 
free boundary. These include the new configurational exponent 71 = ||. 
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The critical behaviour of semi-infinite d-dimensional systems with free surfaces can be 
very rich, with the possibility of "special" or multicritical surface behaviour when the surface 
couplings are sufficiently enhanced [|T],0 . One fundamental model of relevance in this context 
is the semi-infinite two-dimensional n- vector or O(n) model ||. In the limit n — 0, which 
describes self-avoiding walks [|J , the special transition corresponds to the polymer adsorption 
transition [||. 

In this letter we derive the surface critical behaviour of the 0(n) loop model on the 
honeycomb lattice with "ordinary" (O) and "special" (S) boundary conditions on either 
side of a finite strip, i.e. a strip finite in one direction and infinite in the other. Such mixed 
boundary conditions are of current interest for general d, being of relevance to the Casimir 
interaction between two spherical particles in a fluid at the critical point (see, e.g. and 
references therein). In d = 2, with boundary conditions a and b on opposite edges of a strip 
of width N, the free energy per site of a critical system scales as (see |7[ and references 
therein) 

f N ~ f B + f a N~ l + UN' 1 + A afe iV- 2 (1) 

for large N. Here fs is the bulk free energy, f a and /& are surface free energies and A ab is 
the universal Casimir amplitude. This Casimir contribution to the free energy has recently 
been calculated for the O(n) model with various mixed boundary conditions via conformal- 
invariance methods 0. Here we derive this quantity exactly for mixed O-S boundary con- 
ditions, along with the surface scaling dimensions, from the Bethe Ansatz solution of the 
corresponding O(n) loop model || on the honeycomb lattice. 
The partition sum of our 0(n) loop model is defined as 

z = Y.* L y L a a y L b b n p (2) 

where the sum is over all configurations of closed and nonintersecting loops on the lattice 
depicted in figure 1. Here P is the total number of closed loops of fugacity n in a given 
configuration. We take boundary condition a (b) on the left (right) edge of the strip. At 
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n = 0, x is the fugacity of a step in the bulk and y a is the fugacity of a step along the 
left (right) edge. Thus L is the number of steps in the bulk and L a (L b ) is the number of 
steps along the left (right) edge. 

We have found that the equivalent 3-state vertex model can be solved for a number of 
boundary conditions. The possible arrow configurations and their corresponding Boltzmann 
weights are shown in figure 2. Here the phase factors are such that n = s + = —2 cos4A. 
The integrable bulk weights follow in a particular limit of the Izergin-Korepin i?-matrix 
PJlOU, with ts = 2 cos A, or equivalently, with critical bulk coupling [|ll|] 

l/x* = ^2± v / 2^. (3) 



On the other hand, the integrable boundary weights follow |T2[ from appropriate combi- 
nations of the three known reflection or i^-matrices satisfying the boundary version of the 
Yang-Baxter equation [ |13|1 . 



Three inequivalent integrable sets of boundary weights are known to be compatible with 
O(n) symmetry. One set corresponds to 0-0 boundary conditions with boundary 
weights and equivalent critical surface couplings given by 

t a = t b = t B y* a =y* b =x*. (4) 



Another corresponds to S-S boundary conditions 15], with 



cos2A „ „ , N 



where 



l/y* = V±V2^n. (6) 

The third set is a mixture of the above, and corresponds to O-S boundary conditions, with 

cos2A ._, 

t a = t B , h = => y a = x , y b =y . (7) 

cos A 



The self-avoiding walk point at n = occurs at A = 7r/8, where l/x* = y2 + v / 2 
and l/y* = 2 1 / 4 In the lattice model of the polymer adsorption transition 0, the 

self-avoiding walk has energy 



E 



-eh. 



(8) 



where e is a constant and L s is the number of steps along the adsorbing boundary, in this 
case the right hand side of the strip. For the O-S boundary conditions, we thus obtain the 
same critical adsorption temperature 



exp 



{w) = y * /x * = + ^ = L553 • • • 



(9) 



as for the S-S boundary conditions ||15|| , i.e. with adsorbing boundaries on both sides of the 
strip. Recent phenomenological renormalisation transfer matrix calculations on the square 



lattice with one adsorbing and one free boundary are consistent with this finding [16] . 

We have previously solved the corresponding vertex model by means of the co-ordinate 
Bethe ansatz for the 0-0 and an analytic Bethe ansatz for S-S boundary conditions 
H] , [T^|jr5| . Proceeding in a similar manner to [I2|, we find that the eigenvalues of the 



transfer matrix for the O-S boundary conditions are given by 

a _jt smh( Uj + i 3A/2) sinh(w j - i3A/2) 
~~ sinh(w i + i A/2) sinh(w i — i A/2) 

where the Uj follow as roots of the Bethe Ansatz equations 



(10) 



sinhfwj — i A/2) sinh(uj — i 3A/2) 



sinhOj + i A/2) sinh(u j + i 3A/2) 



A' 



sinh(2uj + i A) 
sinh^Wj — i A) 



x 



n 



sinhfwj — u k + i A) sinh(uj + u k + i A) sinhf/Uj — u k — i 2A) sinh(-Uj + -u fc — i 2A) 



^ sinh(-Uj — u k — i A) smlif/u^ + u k — i A) sinh(uj — u k + i 2A) smh(uj + u k + i 2A) 



:n) 



Here N is the width of the strip (e.g., N = 8 in figure 1) and m labels the sectors of the 
transfer matrix, with m = N for the largest eigenvalue A . A more convenient sector label 
is I = N — m. 

The direct calculation of the finite-size corrections to the eigenvalue spectrum follows a 



well-trodden path (see, e.g. JT^| and references therein). In this case the largest eigenvalue 
in a given sector m is characterised by real positive roots with related integers Ij = j, j = 
1, . . . , m. Defining the free energy per site as Jn = A" 1 log A , we find 



f N ~f B + /o-sA- 1 + Ao-sA" 



(12) 



This result is to be compared with equation (1). Dealing with the non-universal terms first, 



r 

h = J_ 



00 sinh( iff — A)x sinh Ax 
xsinh |7rx(2 cosh Ax — 1) 



dx 



(13) 



is the bulk free energy [17, Hj] and 

1 — cos A 



fo-s = 2 lo g(^ 



+ 2 



cos3A 

00 sinh I Ax cosh \ — 2A)x sinh t(vt — 6A)x 



dx 



(14) 



3 xsinh |7rx(2 cosh Ax — 1) 

is the surface free energy for the mixed O-S boundary conditions. We note however, the 
identity 



fo-s — 9 (/oo + /s-s) — fo + /s 



(15) 



where /o-o and /s-s have been derived previously p4]JT5|JT8|] . The individual surface contri- 
butions are thus 



1 / 1 - cos A n 

fo = T lo g(i £t) 

4 v l — cos3A / 

sinh jAx cosh |Ax cosh |(7r — 2A)x sinh 7(71" — 3A)x(2 cosh |Ax — 1) 



+ 2 



xsinh \nx{2 cosh Ax — 1) 



c?x 



(16) 



and 



/ s = i log(— 



1 — cos A 



cos3A 

00 sinh I Ax sinh |Ax cosh |(7r — 2A)x cosh |(7r — 3A)x 



<ix. 



(17) 



x sinh ^7rx(2 cosh Ax — 1) 

In particular, at n = we have A = (2 + \/2) N /(l + v 7 ^) with / B = log(2 + ^2), f o = 
and fs = — log(l + \/2). Here the sign change in fs represents an attraction towards the 
adsorbing boundary. 

The Casimir amplitude appearing in (1) and (0) is given by 



A 



o-s 



24 



(18) 



where ( — 2/ is a lattice-dependent scale factor. The effective central charge is 
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c = 1 - 



12A 2 

7r(7T - 2A) 



(19) 



with c = at n = 0. The mixed boundary scaling index |T9|j20[| follows as 

w ^ -vr + 8A 
*o-s = m( c ~ c ) 



(20) 



8(tt-2A) 

in agreement with the conformal invariance prediction J/J . 

The geometric scaling dimensions defining the surface critical behaviour follow from the 



inverse correlation lengths via 21 



— i 



, A n(X e 
log — ~ 



(21) 



These dimensions govern the geometric correlation 



G e (x-y)~\x-y\- 2X c 



(22) 



between £ nonintersecting self-avoiding walks tied together at their extremities x and y, 
which for surface critical phenomena, are near the boundary of the half-plane |22| . For mixed 
boundary conditions there is a discontinuity at the origin between boundary conditions a 
and b corresponding to the insertion of a boundary operator . 

Here the scaling dimensions are associated with the largest eigenvalue in each sector 
of the transfer matrix. We find 



X e = y£ 2 + Ug-2)£, £ = 1,2, 



(23) 



where rcg = 2tt — 4A. These dimensions are to be compared with the other Bethe ansatz 



results. For the ordinary (O-O) transition [[14 



X e = lg£ 2 + l(g-l)£ 



(24) 



or Xt = hg + i : i in terms of the Kac formula |23| , [22|] . On the other hand, for the special (S-S) 



transition 



X e = lg(£+iy-l(£+l) + 9 [<l " 



(25) 



or Xi = hi + i t 3 p4|-pq]. For the mixed O-S boundary conditions we have Xg = h 



£+1,1 



1. 



At n = (g — |), corresponding to mixed adsorbing and free boundaries in the polymer 
problem, (|23| ) gives 

X e = If - \L (26) 

The first two values are X\ — \ and X 2 = X e = 1. These dimensions define critical exponents 
for polymers in the upper half-plane, with one boundary condition on the positive and the 
other on the negative x-axis, the two geometries (the strip and the half-plane) being related 
via a conformal map (see, e.g. Ref. |19]]). In particular, the number of self-avoiding walks 
which begin near the origin (where the boundary conditions meet in the half-plane) scales 
as L 7l_1 yU L where /i = 1/x* and the universal value 71 = || follows from the usual scaling 
relation 

7i = (2-X 1 -X 1 bulk )z/ (27) 



where v = | and Xj 51 ^ = ^ [ITTJ. The exponent 71 = || is to be compared with the 



exact values for the non-mixed cases, where 71 = || for a non-adsorbing boundary (ordinary 
transition) and 71 = || for an adsorbing boundary (special transition). 
A detailed account of our results is currently in preparation f[8] . 
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FIGURES 



FIG. 1. The open honeycomb lattice. The transfer matrix acts in the vertical direction. 
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FIG. 2. The allowed arrow configurations and corresponding Boltzmann weights for (a) bulk 
and (b) surface vertices. 
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